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Abstract 

In this article we describe the generation of the evanescent waves 
which are present in the rarer medium at total reflection by using a 
mixed-type system, the Ludwig system, which leads naturally to con- 
sider a complex-valued phase. The Ludwig system is derived from the 
Helmholtz equation by using an appropriate modification of the sta- 
tionary phase procedure: the Chester, Friedman and Ursell's method. 
The passage from the illuminated to the shadow region is described 
by means of the ray switching mechanism based on the Stokes phe- 
nomenon applied to the Airy function. Finally, the transport system 
connected to the Ludwig eikonal system is studied in the case of linear 
wavefronts and the existence of the Goos-Hanchen effect is proved. 

AMS Subject Classification: 78A05, 35M10, 34M40 

1 Introduction 

The study of evanescent waves attracts nowadays the attention of many au- 
thors in view of their great relevance in several physical and technological 
problems ^ E] • In this article we want to reconsider what is probably 
the first and more elementary process, which generates evanescent waves: 



the total reflection at the interface between an optically denser and a rarer 
medium. In spite of the fact that this phenomenon has been deeply stud- 
ied, starting from Newton still some open theoretical questions remain. 
These ambiguities emerge with some evidence from what Sommerfeld writes, 
about this subject, in his classical textbook Concerning the propagation 
along the boundary, he writes: "An actual energy flow takes place paral- 
lel to the boundary surface. This seems to contradict both the name "total 
reflection" and our oft-repeated statement that no energy is lost in the pro- 
cess. We must, however, consider the fact that we have always performed 
our calculations for the ideal case of an infinitely wide wave front. For the 
actual, laterally restricted waves, energy can very well pass from the denser 
to the rarer medium or, respectively flow back from the rarer into the denser 
medium at the lateral boundaries of the wave. This is the energy which is 
transported parallel to the boundary surface or, so to say, meanders about 
it". Next, he quotes the experiments of Goos and Hanchen which show 
the penetration of the incident ray into the rarer medium. 
At this point two questions arise: 

(i) Is it possible to find a system which describes globally the propaga- 
tion of rays in the various domains: in the illuminated region (denser 
medium), in the shadow (rarer medium) and also along the boundary? 

(ii) Can this global representation of ray propagation explain the penetra- 
tion of the incident wave into the rarer medium without introducing 
complex-valued transmission angles in a formal way? 

We shall try to answer the previous questions limiting ourselves to the 
scalar representation of the light (which corresponds to linearly TE polarized 
light). In what follows we will describe phenomena which are related to 
the generation of evanescent waves at total reflection, as the plane wave 
propagation along the interface and the Goos-Hanchen effect. Of course, 
this does not exhaust the whole phenomenology; other effects and waves 
which can be present, e.g., surface waves excited at the interface and 
waves exhibiting unusual properties (P* wave) |H1, will not considered here. 

It is well-known that geometrical optics is given by the leading term of an 
asymptotic expansion of a solution of the wave equation, which is usually 
derived by using the stationary phase method. One then introduces the 
concept of ray, which is described by the so-called eikonal equation. In the 
present article we aim at showing that the propagation along the boundary, 
which is a peculiar feature of total reflection, cannot be treated within the 
framework of the standard geometrical optics, but by the use of an extended 
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formulation of the latter based on the Ludwig eikonal system. Indeed, as 
we shall show in Section |21 the standard stationary phase method cannot be 
used to describe the ray propagation along the boundary, when the incident 
ray reaches the critical angle, since two stationary points coalesce at the 
interface between the two media. One is then forced to use a modification 
suggested by Chester, Friedman and Ursell (CFU method) 9 , which leads 
to a non-standard eikonal equation, which was first introduced by Ludwig 
jlUj in a quite general context. 

The paper is organized as follows: In SectionElwe start from the Helmholtz 
equation and derive the Ludwig system by means of the CFU method. In 
Section 121 we discuss in detail the mixed-type (hyperbolic-elliptic) character 
of the Ludwig eikonal system. In Section^the Stokes switching mechanism, 
which allows us to pass from the two-ray illuminated region to the one-ray 
shadow region, comes into play. In Section [21 we study in detail the propa- 
gation of the waves in the various domains, in the case of linear wavefronts. 
We show that the Ludwig eikonal system can, indeed, furnish a global de- 
scription of the phenomena and explain the generation of the evanescent 
waves in the shadow region. Moreover, from the analysis of the Ludwig 
transport system we describe the Goos-Hanchen lateral shift of a reflected 
ray at total reflection jH] , which is related to the penetration of the incident 
ray into the rarer medium. In conclusion, we can give a first, partial positive 
answer to the questions posed earlier in the limit of the scalar treatment of 
the problem. Finally, in the Appendix we illustrate in detail how to derive 
the Stokes switching mechanism by means of the Borel summation method. 

2 Collision of stationary points and the Ludwig 
eikonal system 

Let us start from the Helmholtz equation, which reads 



where k is the free-space wavenumber, and n is the refractive index of a 
homogeneous medium. We look for a solution to of the following form: 



(x E M^). The principal contribution to ip{'x.,k) as /c — > oo, corresponds 
to the stationary point of in the neighborhood of which the exponential 



V V + fc^nV = 0, 
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exp(ifc$) ceases to oscillate rapidly. These stationary points can be obtained 
from the equation 9$(x,/?)/a/3 = (provided that d^^ix, (3)/d(3'^ / 0). Let 
us suppose for the moment that at each point x we have only one stationary 
point of ^, then the asymptotic expansion of -0 as A; — > oo can be written as 
follows: 

oo , 

^(x,fc):.e''=*(-/^o)5:-^, (3) 

m=0 ^ ^ 

where /3o is the unique stationary point of The leading term of expansion 
(jHJ reads 

V^(x,A:):^Ao(x,/3o)ei'^*(^'^°), (4) 



where 

^2$ 



Ao(x,/3o) =^(x,/?o) 



-i/2\ 

exp[i^sgn(9^$/a/3^)^=;3oJ- (5) 

/ /3=/3o 



Substituting the leading term (that hereafter will be written as ^ exp(iA;<I>), 
omitting the subscript zero), into collecting powers of (ifc), and then 
equating to zero their coefficients, we obtain two equations: the eikonal (or 
Hamilton-Jacobi) equation: 

(V$)^ = n^, (6) 

and the transport equation 

V • (^^v^) = 0, (7) 

whose physical meaning is that the current density is conserved. The surface 
$ = const, is called constant-phase surface. 

Now, suppose that two stationary points coalesce. In such a case the 
standard asymptotic method and the corresponding asymptotic expansion 
fail, and a different strategy must be adopted. The procedure suggested by 
Chester, Friedman and Ursell [0] consists in driving the phase function $ 
into a new more convenient form by a suitable change of the integration 
variable /3 <-> ^ (provided $/3/3/3 ^ 0): 

<I>(x,/3)=^/(x) + /(x,e). (8) 

After this change of variable, integral © reads: 

V'(x, k) = e'^"W J 5(x, e) e'*^^^'''^) de. (9) 
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This expression is very similar to integral © with the phase function <I> 
replaced by the function /, and with an additional oscillatory factor exp(iA;u) 
in front. 

Now, suppose that d^/d(3 vanishes at two distinct stationary points 
/3i(x) and /32(x). We want to choose a transformation such that to these 
points there correspond points where df /d^ vanishes, which are symmetric 
with respect to ^ = 0. The result can be achieved by setting 



/(x,0 



1 13 



(10) 



In fact, df/d£, = — .^^ + v(x), and the equation df /d(, = yields ^ = ±\/v. 
Then from equalities (jSJ and (fTI7|) (setting ^ = ^/v and ^ = —\/v), we obtain 
the following relationships: 



$(x,/3i) 
f(x,/32) 



n(x) + 



n(x) 



2„3/2 



(x), 



that yield 



n(x) = i[$(x,/3i) + $(x,/32)], 
x) = i[cI>(x,/3i)-c&(x,/32)]. 



(11a) 
(lib) 

(12a) 
(12b) 



^(x,/32 



In the case Pi = /?2, we have f(x) = and ti(x) = $(x, /3i) 
If equalities (|12|) are satisfied, then the transformation ^ /? is uniformly 
regular and 1-1 near ^ = (see jU])- Prom @ and ()1U() it follows that the 
leading terms in the expression of V'(x, A;), for large k, can be written in 
terms of the Airy function Ai(-) and of its derivative Ai'(-). We have 



ffolxj 



de + /io(x) / 



d^ 



(13) 



where 50 (x) and /io(x) are respectively the first terms of the two formal 
asymptotic series: Em=o fl'm(x)/^'"; and Em=o ^"i(x)/A:™. Then, by using 
the integral representation of the Airy functions we finally obtain: 



V'(x, k) ~ 2^e''="W 



Ai(-A:2/3.) + M^) Ai'( 



/fcV3 



iA;2/3 



-k^/\ 



(14) 



which, substituted into (^, yields a set of four equations. In what follows we 
will focus only on those equations which correspond to the eikonal equation 
(jni), which generate, in a general geometrical context, a complex-valued 
phase and, consequently, evanescent waves. The equations which correspond 
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to the transport equation ^ will be analyzed in Section [51 in the particular 
case of the propagation of linear wavefronts. 

The system corresponding to the eikonal equation reads: 



which is a fully nonlinear, first order partial differential system ruling u and 
V. Hereafter system (|15|) will be referred to as the Ludwig eikonal system 



3 Complex rays and mixed— type systems 

The Ludwig system has been derived in connection with the problem of 
obtaining a uniform asymptotic expansion at a caustic ^01- Let us recall 
that the caustic is the envelope of a family of rays i.e., (i) for every point 
on the caustic it is possible to give a curve of the family which is tangent 
to the caustic on this point; (ii) no curve of the family has a segment in 
common with the caustic jl2| . When total reflection occurs, even at the 
critical angle, at the boundary between the media we do not have a caustic 
in the strict sense of the definition given above. Nevertheless, beyond the 
critical angle, we have propagation along the boundary, and the constant 
phase line of this ray is orthogonal to the boundary. If the latter is the 
horizontal axis, then the wavefront is vertical. It can be reached either by 
varying the constant phase lines orthogonal to the reflected rays, from the 
limiting angle up to the grazing angle, or by varying, in a similar way, the 
constant phase lines orthogonal to the direction of the incident rays. We can 
thus say that the wavefront of the ray running along the boundary is the 
vertical tangent to the cusps generated by the envelope of the wavefronts of 
the reflected rays when it touches, on the horizontal axis, the envelope of the 
wavefronts of the incident rays. Therefore we can regard the boundary as 
the locus of these cusps. We can thus say that on each point of the boundary 
two stationary points coalesce and the standard stationary phase method 
cannot be used. We are then led to the CFU method and, accordingly, to 
describe the geometry of the rays by the use of the Ludwig eikonal system. 
We can then work out the problem through the system (|15|) taking, for the 
refractive index n, a value n > 1 in the denser medium and, for the sake of 
simplicity and without loss of generality, the value n = 1 in the rarer one. 
At the boundary there is a discontinuity in n. 



I Vti|^ + V I Vf 1^ = n' 
Vu-Vv = 0, 



2 



(15a) 
(15b) 



uni. 
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We want now to study the mixed character of the Ludwig system by 
working out the problem in and using Cartesian orthogonal axes, whose 
coordinates will be denoted by x and y. First we write the characteristic 
determinant (see also ^1 ) : 



D 



2ux 


2Uy 


2vvx 


2vv 


dx 


dy 








Vx 


Vy 


Ux 


Uy 








dx 





2 [{uxdy - Uydxf - v{vxdy - VydxY' 



(16) 



Then, the equation of the characteristics for the Ludwig system reads 
f dx^ ^ 



{Uy - VVy) l^-^j + '^{vVxVy - UxUy) + (""x 



0, (17) 



whose discriminant A is 



A = {vVxVy-UxUyf -{Uy-Vvl){ul-Vvl) = v{UyVx-UxVyf = vJ^, (18) 

where J is the Jacobian determinant of the transformation {x,y) <-> {u,v). 
Provided J 7^ 0, and according to the sign of A, we can distinguish among 
the following cases: For v > 0, i.e., A > 0, we are in the hyperbolic case 
(the characteristics are real and distinct); for v = 0, i.e., A = 0, we are in 
the parabolic case (the characteristics are real and coincident); for f < 0, 
i.e., A < 0, we are in the elliptic case (there are no real characteristics). We 
have a system of mixed type [TK] . 

Assuming v > (corresponding to the illuminated region), Eq. ()17|) 
leads to the following characteristic curves: 



dx : dy = {ux ± VvVx) : {uy ± \fv\ 



vv„ 



Then, we can introduce, in the (n, z;)-plane, the curves 
^^^\u,v) = u± = const., 



(19) 



(20) 



which are branches of cubic curves with negative (<!>(+)) and positive ('I''^"-') 
slopes, respectively (see Fig. P). Let us note that at v = 0, these curves 
have cusps with vertical tangent. From (|19j) and (|20j) we have 



dx Ux ± ^/vVx 
dy Uy ± ^/vVy 



(±)' 



(211 
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Elliptic 



U 



Figure 1: The cubic curves (|^ in the (u,w)-plane. 

and, furthermore, it is easy to verify that |V(ii it = n^- Thus, it 

can be concluded that V(u it give the directions of the rays, which 

are orthogonal to the curves = const, (see Fig. 0. These are real rays, 
as real are the characteristic curves of the Ludwig system: for v > 0, this 
system is of hyperbolic type, and the lines = const, are the constant 
phase curves in the {u, w)-plane. 

In the case v = the system is of parabolic type; we are at the boundary 
between the two media and the refractive index presents a discontinuity. 
This situation can be still analyzed as a limiting case: the lines of constant 
phase are the vertical tangents to the cusps generated by the curves ^^^^ = 
const.. The rays are real and orthogonal to the level lines u = const.. 

Finally, let us consider the case v < 0, which corresponds to the shadow 
region and in which the system is elliptic. From ()2U|) we see that the phase 
^^^\u,v) becomes: ^^^\u,v) = u ^ [^{—v)^^"^: i.e., it becomes complex- 
valued. Equation H15b|) informs us that the gradients of u and v are orthog- 
onal; then, where J ^ 0, this amounts to write 



1 




1 


Vy 


Vu\ 




\Vv\ 





(22) 



from which the equalities 

Ux = pvy, (23a) 
Uy = -pvx, (23b) 
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follow, with p{x,y) = |Vn|/|Vi)| = ((l/|Vt;p) -v)^/"^ > (recall that in the 
shadow region n = 1). From (jlSaj) we have that inside the shadow region 
jVup > 1, the equality holding only when |Vf | = 0. Then, it follows that 
p is bounded away from zero, provided that {Vvl is bounded; in this case 
equations form a Beltrami system expressing the fact that the mapping 
(x,y) iu,v) is conformal with respect to the Riemannian metric 

ds'^ = -(dx^ + dy^), (24) 
P 

i.e., that under this mapping an angle a in the (x, y)-plane measured by the 
metric (|24() in taken into an angle a in the {u, u)-plane, measured this time 
in the usual way JSI- Moreover, let us consider a compact subdomain of the 
shadow region defined by < {—vq) ^ (—v) ^ (— ^^i) < +oo; there inside we 
have that p is bounded and bounded away from zero, provided J / and 
|Vn| bounded: i.e., there exists a constant Q such that 

< ^ ^ i + p ^ Q < +00. (25) 

Q p 

In this case system ()23|) is uniformly elliptic, and the transformation (x, y) 
(n, v) maps infinitesimal circles into small ellipses of uniformly bounded ec- 
centricity: i.e., it is a quasi-conformal mapping ^Sl- Moreover, the following 
differential inequality is satisfied: 



^ + Q\J^ (26) 



where J is the Jacobian of the transformation associated to equations H23() . 
Notice that since |Vf ^ = (| Vnp — 1)/(— !»), the gradient of v tends to vanish 
as V increases into the shadow region, so that any oscillation is smoothed 
down by the elliptic character of the Ludwig system for t> < 0. 

Summarizing, we have that for t> > (hyperbolic case) the values of 
the phase are given by = n it for w = (parabolic case) we 

simply have $ = -u; for u < (elliptic case) the complex-valued phase 
= n=Fi|(— w)^''^ naturally emerges; finally, departing from the boundary 
between the two media, that is in any compact subdomain of the shadow 
region, the system becomes uniformly elliptic. 

It is important to stress here that the complex-valued phase emerges as a 
direct consequence of our formulation of the problem, which leads to consider 
the Ludwig system: its elliptic character in the shadow generates a complex- 
valued phase. This situation is deeply different from the approaches of 
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complexification of the real-ray theory, in which the complex values of the 
phase are considered from the beginning and then, by considering all the ray 
intersections with the real space, aim at selecting all the relevant complex 
ray contributions to the solution of the physical problem (see |16| for a 
detailed discussion of pros and cons of this kind of approach) . 

Remark. We have just seen that by fixing u, then v is determined through 
the system ()15() : we have a map n — > f . We can find the map v ^ u as well. 
Both maps can be viewed as Backlund-type transformations which couple 
H15a|) and ()15b|) . For instance, the map v — > u, in the shadow region, can be 
expressed through the following Backlund-type transformation |131 1171 IT^ : 

V. = T^^^?^ ( 1 (. < 0), (27) 



|Vu, 

which couples the two relationships: 

(a) orthogonality of Vu with respect to Vv (see 

(b) the equality: \\7u\^ = 1 - v\\7v\^ (see ((T3aj) ). 



4 Stokes phenomenon and asymptotic behavior 

We can now draw our attention to the amplitude term in solution (|14|) (i.e., 
the term in brackets), and in particular to the Airy function Ai(z), which 
satisfies the Airy equation 

= (28) 

For large \z\, the solutions to (|28j) can be approximated by a linear combi- 
nation of the functions (see the appendix) : 



w 



(a) 



-^/^exp(±|//2)_ (29) 



Evidently w^\z) and w^^\z) are multivalued functions of the complex 
variable z with a branchpoint at 2 = 0; instead, the solutions w{z) of (|28|) 
are entire function of z because the coefficient of w{z) in (|28() . i.e., z, is 
entire. Therefore, when we turn once around the point z = 0, w{z) will 
return to its original value, but and w^^^ will not. It follows that if a 
specific solution 'w{z) of p8|l is approximated at z 7^ 0, by a linear combi- 
nation ciw^^ + C2W^^\ then it cannot be approximated by the same linear 
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combination at ze , that is after one turn around z = 0. The concept of 
approximation here involved must be domain-dependent: different approx- 
imations live in different angular sectors of the complex z-plane. This is 
basically the Stokes phenomenon jl9l I2UI I21j : constructing an uniform ap- 
proximation leads to retain small exponentials though their numerical value 
is dominated by that of large exponentials; across Stokes lines, the multiplier 
of the small exponential changes smoothly but rapidly ^19'. '221 OH], and the 
aim is to establish the relation between multipliers in the different parts of 
the complex domain. These facts are the price to pay when one wishes to 
represent entire functions by multivalued functions. At first sight, it could 
appear strange to represent entire functions by multivalued functions, how- 
ever it is worth noting that the representation by multivalued functions is 
usually the only way in which the wave character of a solution can be made 
clearly explicit; this is the main reason of the pervasiveness of the Stokes 
phenomenon in many different contexts. 

In the case of the approximation of solutions to second order homoge- 
neous differential equations, and in particular of the Airy equation, which is 
our concern here, the asymptotic expansion of a solution is, in general, a lin- 
ear combination of two complex exponentials of the type exp((l)±(z)), {z G 
C). The usual prescription states that Stokes lines, i.e., the lines where 
nonuniformity in the asymptotic expansion occur, are located where the 
phases of the two exponentials are equal, i.e., when Im[(/>-(_(z) — = 0. 

However, not all the equal-phase lines are Stokes lines and, in general, the 
problem is to decide which of these lines are actually Stokes lines. In the Ap- 
pendix, these and other questions related to the Stokes phenomenon are dis- 
cussed extensively; in particular, we show how to locate Stokes lines without 
invoking the equal-phase argument above, and, accordingly, we construct 
the asymptotic expansion of the function Ai(z) in the different angular do- 
mains, determining where and how ray-fields are switched on and off (see 
formulae (|55)) and the related ray structure in Fig. The Stokes phe- 

nomenon plays a central role in our analysis in view of the fact that we must 
connect the illuminated region, i.e., v > 0, with the shadow, i.e., v < 0; this 
connection can be achieved by investigating what happens to the asymptotic 
approximation of the Airy function when we turn around the point v = 0, 
with V regarded as a complex variable. The Ludwig approximation (|14|) 
contains the Airy function Ai{—k'^/^v) (in this approximation we neglect 
the first derivative of the Airy function, and this is legitimate in a region 
close to the boundary between the two media), whose ray structure is shown 
in Fig. [2J3. At this point we apply the results shown with details in the Ap- 
pendix and we note that for v real and positive we have two contributions of 
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Figure 2: (a): Representation of tlic asymptotic expansion of Ai(z). Stokes lines 
(thin lines indicated by S); anti-Stokes lines (tick lines indicated by AS). Dominant 
(tick lines labeled by 'd') and subdominant (thin lines labeled by 's') rays, (b): 
Representation of the asymptotic expansion of Ai(— fc^^'^w) (see also ^B], Fig. 15 
for an analogous picture). 



comparable amplitude, whereas for v real and negative, i.e., in the shadow, 
we have only one contribution with subdominant character. For v G M"*" we 
have (see the definition of w± in (|5^ ): 

1 



2^ 



\-k'^/^v) + -\w+{-k^l^v)] {v G M+). (30) 



At high frequency (i.e., for k +oo), we can approximate w± with 
obtaining for v £ M"^: 



Ai(-A;2/^ 



1 1 



~i7r/4 



+ — 



V 



1/4 



V 



1/4 



{v G 



(31) 

Plugging the r.h.s. of formula (|31l) in the expression of •ip{x,k) given by 
H14() . and neglecting the second term (which is certainly legitimate for small 
values of v, i.e., near the boundary between the two media JOl) we get the 
following approximation: 



1 e^'"/^ 



^ifc(«-|t,3/2) 



+ e" 



-i7r/2 ifc(w+|t)3/2)' 



{v G 



Proceeding in analogous way, we have for v G 

1 



Ai(-fe2/3, 



2^/^ 



i-k'^/^v) {v G M"). 



(32) 
(33) 
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For sufficiently high frequency we can approximate W- with 'w_ , obtaining: 

Again, plugging the r.h.s. of 1)34(1 in the expression of k), given by (|14() . 
we have: 

Note that n + i|(— v)^/^ is the complex-valued phase. 

Remark. The exponentials which are contained in the functions w^^^ {—k'^^^v) 
are of the type exp(ib|fc(— f )^/^) and, consequently, we have Im0-|-(i;) = 
ib|/c|f l^/'^ sin |(phf + vr). In this case we see that care must be taken in the 
application of the prescription lm[0+(t>) — (j)-{v)] = for the localization 
of Stokes lines; in fact, this would lead to identify erroneously three Stokes 
lines at phi; = ±^,7r. The presence of the Stokes line at phi; = vr would 
imply the existence of a second ray in the shadow, which actually does not 
exist. However, as discussed in the Appendix and also shown in Fig. O), the 
line at ph f = TT is not a Stokes line, so that only the wave with subdominant 
character, given in ((35|l. is indeed present in the shadow. 



5 Generation of evanescent waves at total reflec- 
tion 

In this section we obtain the formulae which describe the wave propagation 
when the wavefronts are represented by straight lines in M^, focusing on the 
total reflection at the boundary between the denser and the rarer medium. 
Then we introduce in the system of orthogonal axes (0;X,Y), and denote 
by X and y the corresponding Cartesian coordinates. The interface between 
the two media is represented by the line y = 0. We treat the three regions 
separately: first, the illuminated region lying in the upper half-plane (y > 0) 
with refractive index n > 1; then the shadow region [y < 0) with unitary 
refractive index; finally, the boundary between the two media at y = 0. 

Besides the Ludwig eikonal system (|15() here we will consider also the 
two equations which correspond to the transport equation 0; following 
the procedure outlined in Section [21 that is substituting Eq. (|14() into the 
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Helmholtz equation and collecting powers of fc, we obtain the following 
equations that govern the amplitudes 50 (x) and /io(x) |10j : 

2Vu ■ Vgo + Augo + 2vVv ■ V/iq + vAvh^ + {Vvfho = 0, (36a) 
2V?J • Vfifo + Avgo + 2Vu • V/iq + Au/iq = 0. (36b) 

Also system H36p . which is linear in g^ and /iq, features a mixed-type char- 
acter similar to the one displayed by the eikonal system (jlSjl : i.e., it is 
hyperbolic where v > 0, and elliptic where u < 0. This characteristic will be 
discussed with more details in the next subsections (the transport problem 
can also be treated by means of the Leontovich-Fock parabolic equation; 
the interested reader is referred to [23 |H] ) • 



5.1 The illuminated region: y > Q 

We take u{x) = nxsinOi, where 6i is the incidence angle computed, as usual, 
with respect to the normal to the interface, i.e., the Y-axis, and n is the 
refractive index of the denser medium. The critical angle 9i, at which total 
reflection sets in, is given by sin^^ = 1/n. First let us consider the eikonal 
system ()15p. From the expression of u{x) we have Ux = n sin 9i, and, from 
Eq. (|15b|) it follows that Vx = 0. From Eq. (|15aj) we have: 

y(^^y = n\os^9,, (37) 

from which it follows that in the illuminated region the sign of v is always 
positive. Integrating 1)3 7|) with the constraint that v{x,y) > for y > 0, we 
obtain, for n = constant: 

Iv^/'^ = ny cos Oi {y > 0). (38) 

Then the constant phase lines ^^{u,v) = n it ^v^^"^ = const, are given, in 
the (x,y)-plane, by ^^{x,y) = xs'mOi ztycos^j = const.. Substituting the 
expressions of u{x,y) and v{x,y) in formula (|32|) . we have: 



^ikn{x sin0i—y cos 0i) _|_ ^— in /2 ^ikn{x sin 0i+y cos 9i) 



(39) 

where Gi{x,y,ei) = y?e^'^/'^(2/3ncos 6ii)i/6(5o(x, y)/yi/6), and the sub- 
script I stands for recalling that we are considering the illuminated re- 
gion. Formula factorizes the wavefunction V/ as the product of the 
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amplitude (Gj/Vk) times the exponentials representing the ray propaga- 
tion: the incident ray exp[iA;n(x sin — ycosOi)] travelling along the direc- 
tion (sin^j, — cos^j) with wavenumber kn, k being the wavenumber in the 
vacuum (recall that a time dependence e~''^* is understood and omitted), 
and the reflected ray exp[ifen(a; sin -f-y cos 9i)] travelling along the direction 
(sin ^j, cos 6*4 ). At the grazing angle, i.e., 9i = 7r/2, the wavefunction is given 
by tl^i = {y/2e~^'^^'^Gj /y/k)e^^^^ , which represents a ray travelling along the 
X-axis, and no reflected ray is present, as it should be. 

Remark. There is a very unpleasant ambiguity in the literature concerning 
the sign of the phase-shift e^^'^^'^ associated with the reflected ray. This am- 
biguity is closely connected with an analogous one encountered in connection 
with the Maslov indexes [J^. In formulae (|,39j) the sign of the phase-shift 
(i.e., the minus sign in e~™^'^) plays a relevant role in order to give a correct 
meaning to the physical displacement of the wavefront. Furthermore, this 
sign follows from our analysis of the asymptotic behavior of the Airy func- 
tion (see formula ()88|) of the Appendix), which, in some sense, completes the 
Dingle's representation. Last but not least, this sign is in accord with the 
following Keller's prescription, given in a different but contiguous problem: 
there is a phase-shift of the form e"''^/^ along a path in the direction of 
propagation; the sign must be reversed in the opposite direction pSj . 

Let us now consider the transport system ()36() . From the expression 
of u{x,y) and v{x,y) we have: Au = 0, Vy = (2n^ cos^ 0^/3?/)^/^, Vyy = 
-{{2n^ cos^ 9iY/^ /{2,yfl^), which substituted into ^ yield the following 
first-order system of partial differential equations: 

^ + = 7-50, (40a) 

ay ox by 

dho , _i/3 dgo 1 

-^ + aiy '^ = --^ho, (40b) 
oy ox by 

where aj = ta.n 9i{^n cos 9i)~^/^ and /?/ = tan 9i{^n cos 9i)^/^ . It is con- 
venient to recast system ()40p in vector form by introducing the unknown 
2-vector (column vector) q(x,y) = {go,ho)^\ then, system (|4U|I becomes: 

qj, + Sq^ = C7q, (41) 

where the 2x2 matrices B and C read 
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The characteristic curves associated with system (|41|) can be obtained from 
the characteristic differential equation j27j 

^ = A± = ±^/7^i = ± tan Oi, (43) 
dy 

where X± denotes the two eigenvalues of the matrix B. The characteristic 
curves are then given by: C^{x,y) = x it y tan 9i = const.. Since the 
characteristics are real and distinct, system (|1T|) (and also system (00))) is of 
hyperbolic type. It is important to notice that the amplitude characteristic 
curves = const, are orthogonal to the constant phase lines <I>~ = const., 
while the curves C~ = const, are orthogonal to the lines = const.. 
System ()41() can be transformed in a more convenient form by introducing 
the new unknown column vector p = {p^^\p^'^^)'^ by the relation q = Fp, 
where T denotes the matrix whose columns are given by the real eigenvectors 
^± of the matrix B, corresponding to the eigenvalues \±. The matrix T 
reads: 



tanOi aiy-^l^ } ' ^^^^ 

Therefore, system (|41|) is transformed in the following couple of differential 
equations for the functions p*^^^ and p^^^ : 

p(,i) + A_p« = (A_ = -tan^i), (45a) 

+ ^ 1 (2) ^ ^ (45^) 

Let us now consider Eq. H45a() : recalling that A_ = dxjdy (see Eq. (031)) 
we see that, along the characteristic curves C"*" = const., Eq. H45a|) reduces 
to the following ordinary differential equation: 

d-p^^^ ^ dx dp^^^ 1 f^^^^ 

dy dy dy dx 6y ' 
whose solution is given by: 

p^^\x{y),y) = Ci y^^^'° on the curves C^{x^y) = const.; Ci = const. . 

(47) 

It follows that the solution to Eq. H45a() will be given by Eq. (|47() times a 
function depending only on the characteristic curve C"*" = const.. With anal- 
ogous arguments we find that the solution to Eq. ()45b|) on the characteristic 
curves C~{x,y) = const, is given by: 

p^'^\x{y),y) = C2 y^^^ on the curves C~{x,y) = const.; C2 = const. . 

(48) 
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Then, the general solution to system (|45|) is given by: 



P^^\x,y) = fi{x + yis.nei)y-^'\ (49a) 
= fi{x-yi&nei)y^/\ (49b) 

where //(•) and f^i') are two sufficiently regular arbitrary functions. 

We can now return back to the vector q and to the functions go and by 

means of the relation q = Fp: we have 



go{x,y) 
ho{x,y) 



^ fiix + ytan^j) + ^{{x -ytanOi) ,(50a) 



,1/6 



% fi{x + ytan^i) + ^/a7/2^(x - y tan( 



(50b) 



Finally, by inserting Eq. H50a|) into Eq. p9() we obtain the expression of 
the wavefunction in the illuminated region (y > 0): 



ipi{x,y,9i;k) 



^5/ f({x + y tan 9i) + ^/aj fi{x -y tan 9i 



^ikn{x sinSi— y cos 9i) _j_ g— i7r/2 gifcn(a; sin^i+j/ cos S^) 



,(51) 



with ip'j = const., and // (•) and (•) are (at this stage) arbitrary functions. 

The field ipj is the product of two factor: the first one, which arises from 
the transport system ()36|1. describes the distribution of the amplitude on the 
lines x±y tan^j, while the second one, which arises form the eikonal system 
(jlSf) , represents the evolution of the constant phase lines x sin 9i^y cos 9i = 
const.. Obviously, in the wave propagation in a dielectric, homogeneous and 
isotropic medium, like the one being considered here, the constant phase 
lines must be orthogonal to the constant amplitude lines; then from ()5ip we 
are led to consider 



i^i{x,y,9i;k) 



+ e 



^ifi {x + ytaxi9i)e 



ikn{x sin 8i—y cos 8i) 



-i7r/2 



aifiix -yta.n9i)e 



ikn{x sin 8 i+y cos 6i ) 



(52) 



Formula (|52() provides the wavefield in the illuminated region as the sum of 
two terms: an incident wave propagating along the direction (sin 9i, — cos 9i) 
whose amplitude on its linear wavefronts is described by the function //, 
and a reflected wave propagating in direction {sin 9i, cos 9i) with amplitude 
distribution described by /2 on its linear wavefronts. 
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5.2 The shadow region: y < 

In agreement with the expression of u{x) taken in the illuminated region, we 
now set: u{x) = s'mOtx (recall that the refractive index in the region y < 
is unitary). But, now, the incident and the refracted rays invert their role; 
therefore we can still use the Snell formula but putting 1/n in place of n: 
i.e., writing: sin^j/sin^j = 1/n. Therefore, we have sin^^ = nsin^j, and, 
accordingly, u{x) = nsinOiX. From Eq. ()15b|) it follows that Vx = 0, and 
from Eq. ()15a|) we have: 



p-Y = l-n^sin^9. = l-'^, (53) 



from which it follows that in the region y < the sign of v is negative for 
9i > 9^, i.e., when total reflection occurs. In what follows we will assume 
9i > 9i. Integrating (|53|) with the constraint that v{x,y) < for y < 0, we 
obtain, for n = constant: 



(-i;)3/2 = sin2 9i - 1 y {y < 0] 9i > Og, V < 0) . (54) 



Substituting the values of u and v in formula we obtain: 



where Gs{x,y,9i) = ^(|^/n^^l?9~^)-^/<^{go{x,y)/{-y)^/^), and the 
subscript S stands for recalling that we are studying the propagation in the 
shadow region. Formula ()55() factorizes the wavefunction as a product of 
an amplitude {Gg/Vk) times an evanescent wave which propagates along the 
positive direction of the horizontal axis, and which vanishes exponentially 
along the vertical axis. 

Let us now consider the system (|36p . From the expressions of u{x,y) 
and v{x,y) we obtain: Au = 0, Vy = -[2(n^ sin^ 6*1 - l)]-^/^(3y)"-'^/^, Vyy = 
[2{ri^ siii^ 9i - l)]l/3(3y)"'^/^ which, inserted into give the following 

system of partial differential equations for the functions go{x, y) and /io(x, y) 
{y < 0): 

dgo . r, 1/3 1 

+ = 7-50, (56a) 

oy ox by 

9ho . -1/3 9go 1 

—- + asy ' = -7-^0, (56b) 

oy ox oy 
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whereas = (2/3) (n^ sin^ 6* — 1)^2/3^ g^^^^. ^.nd Ps = -{3/2)^/^{n'^ sin^ Oi- 
l)~^/^nsin^j. Although the structure of this system is very similar to that 
of system ()4U[) obtained in the illuminated region, the two systems are pro- 
foundly different. In fact, the characteristic differential equation for system 
H56jl reads: 



dx\' 



v? siv? ( 



sm 



11? sin^ Oi 



< 



> 



sm 



(57) 

so that system does not have real characteristic curves: it is elliptic. 
Let us introduce the functions goix^y) and /io(x,y), defined by: 



ho{x,y) 



{-y)'/%o{x,y). 



(58a) 
(58b) 



Substituting Eqs. (|58|) into system (|56|) and after suitable differentiation, 
we obtian the following couple of differential equations: 



dy"^ 



dx"^ 



0, 



dy"^ 



asPs 



dx'^ 



0, 



(59) 



which are Laplace-type equations {asPs < 0), whose general solution can 
be written as the sum of two arbitrary analytic functions of the complex 
variables {x + '\\/—asPsy) and {x — W—asPs u), respectively pS'. Thus, 
from the solutions of equations (|59() and taking into account Eqs. H58|) and 
system we obtain, for y < 0: 

go{x,y) = i-y)^/^ \ffix + W-asPs y) + fi{x - \^/-asPs y)] ,(60a) 



as , 



/io(x,y)=iJ-^(-y)-V6 



/f (x + i^/-asPs y) - fl{x - W-asPs y) 



(60b) 



By inserting Eq. (jbUaj) into Eq. H55|) we finally have the expression of the 
wavefunction for y < and 0i > Of. 



r 



ips{x,y,Oi;k) ~ /f (x + i\/-asPs y) + fiix - iy'-asPs y) 



X e 



iknxsindi pfc\/ sin^ 9i — ly 



(61) 



with = const., and /f (•) and /2'(') being arbitrary analytic functions. 



19 



5.3 The boundary between the two media: y = 0, and the 
Goos— Hanchen effect 

At the interface y = the refractive index is not defined since it presents 
a discontinuity, and, therefore, the eikonal and transport systems can be 
analyzed only in the limit of y — > 0^. From the expressions of v{y) in the 
illuminated and shadow regions, we see that v{y) — > for y ^ 0. Then, in 
the limit sense, the l.h.s. of Eq. (|15a|) reduces to jVup, and we reobtain, 
in a certain sense, an eikonal equation in standard form, except that at the 
r.h.s. of (|15a)) is not well-defined. 

Since the wavefield and its normal derivative must be continuous, the 
following boundary conditions at the interface y = must be satisfied: 

ipi{x,y)\y=o+ = i;s{x,y)\y=o-, (62a) 



along with appropriate radiation/decay conditions at infinity. 

The continuity constraint ()62a|) yields the following relationship among the 

functions //, f^, ff and f^: 

i^s [/f (^) + /2''(^)] = V'? \v%fl{x) + e--/2y^/|(x)l . (63) 



From the boundary condition (|62b|) . and taking into account the relation 



H63() . we obtain the following relation between the functions // and 




cos Ui — IV sm Uj — sm 



cos tli + i ysin^ 9i — sin^ ^ 



fli^)- (64) 



The term in brackets is of the form z/z* , i.e., it is a phase factor; then ()64|) 
can be written as 

/^(x)=e-/^/^e^^(^«)//(x), (65) 



aj 

where 



5{9i) = — 2arctan — ^- — ^' — — = — 2 arctan 4 /— — (66) 

cos Oi y as PS 

By inserting ()65() into H52|). we obtain the following expression of the field 
in the illuminated region (y > 0): 

'4ji{x,y,9i;k) ~ ^ //(x + ytan0i)ei'="(^'^''^^»-2^^°^^») 
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+ e^^(^'V/(2; - 2/tan0i)e''="(^''^'°^'+^™'^^')j , (67) 

= const.. Finally, condition (|62b|) and formula H65j) yield the following 
expressions for the functions //(x), /f (x): 

^^/f(x) = V/V^[l + itan(5/2)]//(x) + C, (68a) 
V;^/|(a;) = V'/V^ [1 + itan(V2)] cos 5 //(x) - C, (68b) 

C being an arbitrary constant. Plugging equations (jKH]) into (|HT|) . we finally 
have the expression of the field in the shadow region (y < 0): 

ijs{x,y,Oi;k) ~ ^[l + itan(V2)] 

f({x + iy/-asPs y) + cos (5 fi{x - \yj-as(is y) 

^ gifcnxsinei ^kyj 'n? sin^ 8i-ly ^-gg-j 

The reflection coefficient e'^^^'^ in (|67jl describes the phase change that the 
reflected ray undergoes at total reflection, and coincides with the phase 
change given by the Maxwell's theory j29| I.SOj for a linearly polarized TE 
wave (also called s polarization). It is worth noting that this phase-shift is 
intimately related to the mixed hyperbolic-elliptic character of the Ludwig's 
transport system ()36() since, as shown by the formula of 6{6i) in ()66() . it 
depends on the ratio between the characteristic determinants of system (|36|) 
in the illuminated and in the shadow regions. 

The role of e'"^^^*^ can be easily illustrated. The wavefront associated with 
the reflected ray is delayed by the factor e^^ (notice that (5 ^ 0), so it appears 
as if it would have been reflected by some virtual surface located a small 
distance within the rarer medium and not by the actual physical interface 
(see Fig. ISJ. We thus find a behavior already conjectured by Newton [1]: 
at the condition of total internal reflection at the boundary between an 
optically denser and an optically rarer medium, a beam of light penetrates 
some distance into the rarer medium, and then re-emerges in the optically 
denser medium. Consequently, the actual reflected (i.e., totally reflected) 
ray is laterally shifted with respect to a ray geometrically reflected (see Fig. 
121). This type of behavior has been experimentally verified by Goos and 
Hanchen indeed jH], and one can thus speak of Goos-Hanchen effect |6l I31j. 
The magnitude of this lateral shift, which is found experimentally to be of 
the order of a few wavelenghts, can be evaluated as follows [H^ OHl OH 155] . 
Since (5 is a function of the incidence angle 9i, it can be thought of as a 
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Y Virtual line 



Figure 3: Path of a light ray at total reflection. The incident ray is not geometri- 
cally reflected at the interface y = 0, but at the virtual line located at y = Y into 
the rarer medium, so that it emerges shifted by D in the denser medium. 



function of the x-component kx = kn sin 6i of the wavenumber vector; then, 
expanding 6{kx) in a Taylor series about kx = 0, and truncating at the hnear 
term, we have: 



S{kx)=So + Xkx + 0{k',), X 



^ A 



1 d5 



27rn cos Oj dOj 



^ 0, 



(70) 



where A is free-space wavelength. Then, substituting (jTUf) into formula ()67|) 
we have 



tpj{x,y,9i;k) ~ 



f({x + ytan0i)e 



ikn{x sin 9i—y cos 6i) 



+ e'^'^flix- y tan 9i)e 



ikn[{x+X) sin 9i+y cos 9i] 



(71) 



that is, the reflected ray appears shifted by \X\ in the positive x direction 
(recall that X < 0). From the geometry shown in Fig. |3J the measured 
Goos-Hanchen lateral ray-shift D can be easily computed; we have: 



D = X cos ( 



A d5 
2Tm dOi ' 



(72) 



which is identical to the classical expression obtained by Artmann in [32 
(see also jlS]). 
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Finally, by inspection of (|51() and in view of formula H55j) relating f[{x) 
and flix), we can see that, if the function // is constant, then the factor 
in the first square brackets in H51() becomes merely a complex constant, and 
the phase-shift between the incident and the reflected ray reduces to e~^'^/^: 
i.e., 5{9i) does not play any role. Then the Goos-Hanchen shift due to 
5{9i) appears only when the amplitude on the wavefront is not constant, in 
particular when the incident wave is a light beam with limited transverse 
extent (see also [311 1^ IMl H] for a discussion of the Goos-Hanchen effect 
in the case of incident beams with finite transverse width). 



Appendix: Borel summation and the Stokes switch- 
ing mechanism 

The main purpose of this appendix is to show how to locate Stokes and anti- 
Stokes lines without invoking any physical, though reasonable, argument, 
but working out the problem exclusively by the use of mathematical tools. 
We split this appendix into three parts. The first part is devoted to the Borel 
summation which, after Dingle ^Q^, is the main mathematical tool used for 
investigating the Stokes phenomenon. The second and the third parts are 
devoted to the discussion of the Stokes and anti-Stokes lines, respectively, 
(a) Borel summation. Consider the formal series: X]jS=o(~-'-)"("''/-^")- 
Following the Borel method, we formally write: 



n=0 ^ 



+ CXD 



E 



.n=0 ^ ^ 



dt. (73) 



We can give a meaning to the integral at the r.h.s. of ()73() in view of the 
following considerations: the series '}2'n=o^~^ I ■^)^ converges, in the domain 
\t\ < X {x > 0), to a function 1/(1 + (t/x)), which can be continued analyti- 
cally outside the region \t\ < x to the whole semiaxis [0, +00), and which is 
dominated by the exponential e~*. Therefore the integral 

f+00 -t 

converges, and represents the Borel sum of the formal series at the l.h.s. 
of (|73|) . Next, we can perform the analytic continuation of Ao(x) from 
to C. The integrand of the function Aq(z) presents a pole at t = —z, and 
therefore, when one turn is made around the origin, we get: 

Ao(ze^^^') = Ao(z) ^ 27^ize^ (75) 
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where the factor e^^'^' is for indicating that one turn is made around the ori- 
gin in counterclockwise (+) or clockwise (— ) sense, respectively. Proceeding 
in an analogous way, we can perform a summation, in the sense of Borel, 
not of the whole series but only of the late terms. We have: 



E 



dt 



{-xy 



-oo p g-t 



1 + {t/x) 



We are thus led to write: 

oo 

m=n 

where A„(x) is given by 



ra\ 



X" 



Ari(x), 



dt. 

(76) 

(77) 



1 



An (a;) = -r 



-oo g-t 



dt. 



n\ Jq 1 + (t/x) 
Next we obtain a formula, analogous to (|75() . of the following form: 

„±inn (n+l) z 

A„(ze±2-)=A„(z)T27ri- 



(78) 



(79) 



(b) Stokes switching mechanism. Let us come back to the Airy equation 
(|28jl : applying the WKB ansatz, that is assuming that 'w{z) is of the type: 



w{z) = h{z) exp / 0^1'^ dQ ) , 
we obtain the following approximations: 



(a) 



V^exp ± / C'/^dC 



-1/4 

z ' exp 



(80) 



(81) 



The complete asymptotic expansions of two solutions to the Airy equation 
are given by: 



w± = z 



-1/4 



exp (±§2^/^) W±{z) {z e M+), 



(82) 



where the corrective terms VF±(2:) are represented by the formal asymptotic 
series ^[]°(ibl)'"Wm- Using standard notations we have: 

A\{z) = -^w.{z) B{{z) = -^w+{z) (phz = 0), (83) 
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and we now focus on Ai{z). Then we have (see |2U|'): 



n-l 



W.{z) = ^{-irwUz) = ^{-irWm{z) + ^{-irW^{z). (84) 

m=0 m=0 m=n 

Dingle [20] has been able to derive the expression of the late terms Wm{z) 
in the sum at the r.h.s. of (|84() . One has: 

Wm{z) = ^^^^ J}~^^'^^^ {m-s- 1)! Ws{z) (m » 1), (85) 

where J^{z) = \z^l'^ ■ Substituting (jHSJ in X]TO=n(~l)™^m, and using for- 
mula (|77j) . we obtain: 

oo ^ 

^ (-ly^VFn. = j;(-^)^(n-.-l)! VF, A„_,_i(^) (Iph^l < vr) 

(86) 

Let us now explore if there exist and where the Stokes lines associated with 
W-{z) are located (see Fig. |5^). Recalling that the integrand in k^iz) 
presents a pole at t = — which, according to formula (jTH]). generates a 
discontinuity of the form =F2i7r(e^''^"z*^"'"^^'*e^) /n!, we can say that the unique 
Stokes line, associated with VF_(z), lies at ph.F = tt (see formula (|HH|)). i.e., 
at phz = Ivr, recalling that T = Next we apply directly formula 

(|79jl to A„_5_i(.7^); this amounts to take into account a term of the form 
1YK{—T^^"'~^^e^l (n — s — 1)!. Therefore, the asymptotic expansion of Ai(2;), 
which is given by the first equality in 1)83^ when phz varies in the range 
^ phz < Ivr, changes, crossing the Stokes line at phz = |7r, and we have: 



Ai(z) = 

^ ^ 2^ 



m=0 



= (u;_ + it/;+) . (87) 

At this point one can easily verify two of the Dingle's rules j2Uj : 

(i) at phz = Ivr (i.e., ph.F = vr) all the late terms of W-i^z^ are homoge- 
neous in phase and all of the same sign; 

(ii) crossing the Stokes ray at phz = |vr, ^Y^=rX~^y^^rn generates a 
discontinuity, according to formula ((7^ . which is, on the ray, | out of 
phase with the series, and proportional to its associated function. 
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Representation 1)8 7|) contains the function w^, and, therefore the asymp- 
totic series W+. Then, in order to obtain its range of vahdity, we have to 
find the Stokes Unes associated with W+. From the definition of and 
by means of arguments similar to those used above for W-, we have that 
the Stokes Unes associated with are located at the phases for which 
ph {—J-) = vr, that is ph (J^) it vr = vr; this means that the Stokes lines of 
Ty+ are located at phz = 0, |7r. It follows that representation (|87j) holds in 
the range |7r < phz < |7r. 

Finally, since both W± have no Stokes lines in the range |7r < phz < 27r, 
the representation of Ai(z) in H83|) can be extended to this phase range. 

Summarizing, we have: 

^.^^^ _ 1 ^( w^{z) (0 < phz < Ivr) U (Ivr < phz < 27r), 

2^/^T I w^{z) + iw+{z) |vr<phz<|vr. 

(88) 

It is worth remarking that the Stokes lines are associated with the asymp- 
totic expansion of a particular solution to the differential equation ()28() and 
not with the differential equation itself. For instance, to the expansion w^{z) 
is associated the Stokes line at phz = |vr, while the lines at phz = 0, |vr are 
connected to the expansion w+{z). Then, although the asymptotic expan- 
sion of Ai(z) is in general a linear combination of w+{z) and W-{z), Ai{z) 
does not have all the Stokes lines associated with w±{z); in fact, the Stokes 
line at phz = 0, which pertains to w+{z), does not play any role in the 
uniform asymptotic expansion of Ai(z) (see Fig. |2^). 

(c) Anti— Stokes lines. The modulus of 'w±{z) is proportional to the func- 
tion exp [lb cos (Iphz)] . The anti-Stokes lines are the lines delimiting the 
phase range in which one asymptotic series is exponentially dominant on 
the other, while on the anti-Stokes lines the two contributions are com- 
parable. Then the anti-Stokes lines are located at those phases such that 
cos (Iphz) = 0, that is phz = ±f ,vr. In the phase range phz G (— f ) f) 
the asymptotic expansion of Ai(2:) contains only the term w-{z), which has 
a subdominant character since cos (Iphz) > 0. In the ranges |vr) and 
(— ^, — |vr) the Airy function Ai(z) is still represented by the sole term w-{z) 
which, however, acquires a dominant character since cos (|phz) is negative. 
We can thus infer that crossing the anti-Stokes lines at phz = it^, w^{z) 
varies from a subdominant to a dominant character. As we have shown in 
part (b) of this appendix, at phz = if^r the term w+{z) emerges and Ai{z) 
is represented by the combination {w-{z) + iw-^{z)). The term W-^{z) is 
proportional to cos (Iphz) and has a subdominant character in the ranges 
(|vr,vr) U ((— |vr, — vr) since, when phz varies in these intervals, cos (|phz) 
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is negative. For the same reason, w-{z) (proportional to — cos (|phz)) is 
dominant. Let us only note that at phz = vr wc have two contributions 
of comparable amplitude, whereas at phz = we have only one term with 
subdominant character. 
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